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Abstract — A successfully designed bike should possess safety and comfort for the riders. A safe bike means that its 
structure must be strong enough to prevent from damage due to various external loads and a comfortable bike means that its 
suspension systems must be excellent enough to reduce the transmissibility of disturbance coming from the uneven roads to 
the rider. In order to achieve the above goals, various methods have been presented; however, most of them assumed that 
each part of a bike is a “rigid body” except the helical (coil) springs. For the last reason, this paper tries to use more 
versatile finite element method (FEM) to perform the static and free vibration analysis of a bike. It is believed that a finite 
element model with all parts of a bike replaced by the “elastic” elements or lumped masses should be more realistic. In this 
paper, the entire bike structure is modeled by using three kinds of beam elements: pinned-pinned (P-P), pinned-clamped (P- 
C) and clamped-clamped (C-C) beam elements. Among the main parts of a bike structure, the mainframe and rim are 
modeled by the C-C elements, the elastic effect of each tire is modeled by using a P-C element, and each spoke or each 
“ spring-damper unit” is modeled by a P-P element. The key point of this paper is to study the influence of some pertinent 
parameters on the lowest several natural frequencies and mode shapes of the bike. It is found that the radius of the hub 
(disks), the pretension of each spoke, the mass of various attachments or rider, and the riding gesture of a rider have 
significant influence on the free vibration characteristics, the static deformations and internal forces (and moments) of the 
pertinent structural members of a bike. Because the mass of a rider is much greater than that of the bike structure itself, the 
static and dynamic characteristics of a bike with and without a rider on it must be studied, separately. 

Keywords — Rider, Bike, Suspension System, Finite Element Method, Elastic Element, Lumped Mass, Main Frame, Rim, 
Tire, Spoke, Natural Frequency, Mode Shape 


I. INTRODUCTION 

In Ref. [1], Champoux et al. have indicated that “the more manufactures can learn and understand about the dynamic 
response of their products, the more they will be able to benefit both current and potential riders”. It is the last reason, some 
researchers have devoted themselves to the study of vibration characteristics of bikes [1-4], Besides, under the assumption 
that each part of the bike is a “rigid body” except the helical (coil) springs, some researchers paid their attentions to the 
design of rear suspension system of mountain bikes to improve the riding performance and rider comfort [4-8]. Since the 
conventional finite element method (FEM) with the entire bike structure replaced by a number of “elastic” members and 
lumped masses is more able to model a bike “accurately” and “realistically”, and it can be used to study both the dynamic 
and static characteristics of bikes appearing in the foregoing literature [1-8], the objective of this paper is to continue the 
static and dynamic analyses of mountain bikes of Ref. [4] by using the FEM. 

In general, an entire bike is composed of several sub-systems, such as the power transmission system, speed adjusting 
system, braking system, turning system, main frame and wheels. Among various parts of the entire bike, only those affecting 
its stiffness matrix are modeled by finite elements and those contributing to the overall mass matrix only are considered as 
the lumped masses attached to the associated nodes. Based on the last concept, the particulars for the finite element model of 
the bike studied in this paper are stated as follows: (i) The head tube, the front fork and the main frame (including top tube, 
down tube, seat post, seat stay and chain stay) are modeled by using the two-node clamped-clamped (C-C) beam elements 
[9], and a few of the last C-C beam elements are replaced by the pinned-clamped (P-C) ones [10] if one of two nodes of a C- 
C beam element is pinned, (ii) The front and rear rims are similar to the circular rings, thus, each on them are modeled by 
using a number of two-node C-C straight beam elements [11]; for convenience, the total number of beam elements for each 
rim is taken to be the same as that of the spokes on each rim. (iii) Since the stiffness of front or rear hub is much greater than 
the stiffness of each of the attaching spokes, thus, either front or rear hub is assumed to be a rigid body and each spoke 
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connecting a hub and the associated rim is modeled by a pinned -pinned (P-P) beam element, (iv) The total mass of each tire 
is uniformly distributed along the circumference of the associated rim and considered as part of the rim mass per unit length 
(i.e., the effective mass density of the rim is determined by p r = p nm + p tirc with p rim denoting the mass density of the rim 
material itself and p tire denoting the equivalent mass density of the tire with respect to the volume of the rim). The elastic 
effect of each tire is modeled by using a P-C beam element with its axial stiffness k determined by k = F/8 , where F denotes 
the vertical load on a wheel and S is the vertical deflection of the tire, (v) The “spring-damper unit” in the front or rear 
suspension mechanism is modeled by using a P-P beam element with its axial stiffness k determined by the similar way like 
that of the tire. 

For the dynamic analysis, the main objective of this paper is to investigate the influence of the next parameters on the free 
vibration characteristics (such as natural frequencies and mode shapes) of the bike: (i) the masses of the attachments; (ii) the 
radius of hub (disks); (iii) the pretension of each spoke; (iv) the mass (and riding gesture) of the rider. For the static analysis, 
the influence of the last parameters on the deformations and internal forces (and bending moments) of any structural 
members may be studied. However, to save space, only the influence of the mass (and riding gesture) of the rider on those of 
some pertinent members is studied. 


II. Formulations of the Problem 

The equation of motion for the free vibrations of an un-damped structural system is to take the form, [/«] {«} + [/:]{»} = 0 , 
where [m] is the overall mass matrix, [k] is the overall stiffness matrix, {u} is the overall displacement vector and {m} is the 
associated acceleration vector. Thus, the information required for constructing the matrices [k] and [in] are presented in this 
section. 


2.1 Stiffness and mass matrices for the three kinds of beam elements 

The three kinds of beam elements adopted in this paper are shown in Fig. 1 . Each element has two nodes represented by ® 
and ©, respectively. Fig. 1(a) shows the pinned-pinned (P-P) beam element, there are two degrees of freedom (dof’s) at each 
node. The force-displacement relationship for this P-P beam element is given by [9] 
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In the above expressions, {Si^and {Misrepresent the node force vector and node displacement vector of the P-P beam 
element, respectively, and [k] PP is the corresponding stiffness matrix with its coefficients given by Eq. (A.l) in Appendix A. 
Furthermore, the symbols p , E, A and t appearing in Fig. 1(a) represent mass density. Young’s modulus, cross-sectional 
area and length of the beam element, respectively. 
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Fig. 1(b) shows the pinned-clamped (P-C) beam element, there are two dof’s at node © and three dof’s at node ©. The force- 
displacement relationship for this P-C beam element is given by [10] 

{>S}p C — [k] PC {u} PC (3) 

where 
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{u} pc 

= 1© 

u 2 u 

U 4 

u 5 

r 

(4b) 



k a 

k \ 2 

k\3 

k l4 

k\5 




k 21 

b 

^22 

b 

/v 23 

^24 

^25 


Iklpc 

= 

k n 

b 

^32 

b 

^33 

k 34 

k$5 

(4c) 



k 41 

b 

/A 42 

k 43 

k 44 

k 45 




^51 

k^2 

ks3 

k$4 

L/x 

c/i 





where { S } PC , { u } PC and \k | (1( . represent the force vector, displacement vector and stiffness matrix of the P-C beam 
element, respectively, and the coefficients for [/r] PC is given by Eq. (A. 3). In addition to the symbols, p , E, A and i , have 
been defined previously, the other symbol appearing in Fig. 1(b), /, represents the moment of inertia of the cross-sectional 
area A. 

Fig. 1(c) shows the clamped-clamped (C-C) beam element, there are three dof’s at each node. The force-displacement 
relationship for this C-C beam element is given by [9] 


where 
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Where {5} cc , { u } cc and [ k | cc represent the force vector, displacement vector and stiffness matrix of the C-C beam 
element, respectively, and the coefficients of [ k ] cc is given by Eq. (A. 5). 
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In addition to the above-mentioned stiffness matrices defined by Eqs. (A.l), (A.3) and (A.5), free vibration analysis of a 
structural system also requires the mass matrix for each of the constituent members. The mass matrices, [m ] PP , [m] PC and 
[m ] cc , for the P-P, P-C and C-C beam elements are given by Eqs. (A. 2), (A.4) and ( A.6 ) in Appendix A, respectively. 



Fig. 1 The three kinds of two-node beam elements adopted in this paper: ( a ) P-P, ( b ) P-C and ( c ) 

C-C BEAM ELEMENTS 

2.2 Transformation matrices for the three kinds of beam elements 

The element stiffness matrix [k\ q and mass matrix [m] q (q = PP, PC or CC) introduced in the last subsection are obtained with 
respect to the local coordinate system xy. In the conventional FEM, the overall stiffness matrix [A] and mass matrix [M] 
are obtained from the element stiffness matrix [k ] and mass matrix [m] q for each of the structural members composed of 
the entire structure, by using the numerical assembly technique. Where \_k\ q and [ m ] q are the property matrices with respect 
to the global coordinate system xy . In other words, each element stiffness matrix [k\ q and mass matrix [m] q must be 
transformed into [ /c | and [ m ] q by using the following formulas 
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[k\=m T q [k\ q m q 
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then they may be used for assembly. 

In Eq. (7), the symbol [A] denotes the transformation matrix for the <7- type beam element. For the P-P beam element as 
shown in Fig. 2, its node displacements with respect the local xy coordinate system, u t (i=l— 4), and the corresponding 
ones with respect to the global xy coordinate system, (/ = 1-4), have the following relationship 


{ll} pp {if } PP 
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where {;<} = [»] u 2 m 3 u 4 f and {w} = [h 1 u 2 w 3 » 4 ] 7 represent the displacement vectors with respect to the local xy 
and global xy coordinate systems, respectively, while [X\ PP is the transformation matrix of the P-P beam element given by 
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Fig. 2 The node displacements of a P-P beam element with respect to the local xy coordinate 
SYSTEM, U‘ ( i = 1 - 4 ), AND THE CORRESPONDING ONES WITH RESPECT TO THE GLOBAL xy COORDINATE 

SYSTEM u f (t =1 — 4). 


where 


k = cos 9, /t = sin 0 


(10a,b) 


In the last expressions, 9 is the angle between positive x-axis and positive x -axis as one may see from Fig. 2. In practice, the 
values of v and // are determined by 
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where (AJ , v, ) and (x 2 , y 9 ) are the global coordinates of node © and node © of the beam element as shown in Fig. 2, 
respectively. 

Similarly, the transformation matrix for the P-C beam element is given by [10] 
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C beam element is given by [9] 
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2.3 Geometric stiffness matrices for various beam elements 

In order to investigate the effect of pretension T s in each of the spokes on the free vibration characteristics of a bike, the 
geometric stiffness matrix | k G ] (/ for the relevant beam elements are presented in this subsection [12], 

For a spoke modeled by using the P-P beam element, the pretension T s in it will increase its axial stiffness. Thus, its 
geometric stiffness [k a ] pp is given by 
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where t s represents the length of each spoke. 

For a spoke modeled by using the P-C beam element, the pretension T s in it will increase its stiffness associated with both 
the axial dof’s and transverse dof’s. Thus, its geometric stiffness \k G ] PC is given by 
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From Fig. 1(b) one sees that, for a P-C beam element, its 1st and 3rd dof’s are in the axial direction, and its 2nd and 4th dof’s 
are in the transverse direction, thus, in Eq. (15), the matrix coefficients k Gmn with m.n = 1 or 3 represent the contribution of 
pretension T s to the stiffness of axial dof’s, and those with m,n = 2 or 4 represent the contribution of pretension T s to the 
stiffness of transverse dof’s of the P-C beam elements. Furthermore, since the pretension does not affect the stiffness of 
rotational dof’s with displacement u 5 (cf. Fig. 1(b)), the matrix coefficients k G m5 and k G5n with m,n = 1-5 are equal to 
zero in Eq. (15). 

When there exists pretension T s in each spoke, the front or rear rim will subjected to the compressive force T r . Since the rim 
is modeled by a number of C-C beam elements, the geometric matrix [k G ] cc for each of the rim elements is given by 
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where l r represents the length of each rim element. It is noted that, in Eq. (16), all diagonal coefficients are “negative” and 
all off-diagonal ones are “positive”. This result is opposite to Eqs. (14) and (15), because T s in Eqs. (14) and (15) is a tensile 
force and T r in Eq. (16) is a compressive force. Furthermore, for a C-C beam element as shown in Fig. 1(c), its 3rd 
displacement « 3 and 6th displacement u 6 belong to the rotational dof’s, the associated matrix coefficients in Eq. (16) are 
equal to zero, i.e., k c m3 — k G m6 — 0 for m = 1—6 and k G 3n — k G6n =0 for n = 1 — 6 . 

2.4 Rim force induced by pretension of spokes 

For simplicity, in order to determine the rim force T r (appearing in Eq. (16)) induced by the pretension of all spokes in a rim, 
it is assumed that the tensile force T s in each spoke directs to the center H of the hub (cf. Fig. 3). If the total number of 
spokes is N, and the average radius of the rim is r , then the average central force per unit length of the rim is given by 


P = 


N T 

s s 

2 n r r 


( 17 ) 


For the free-body diagram of the half rim shown in Fig. 3, the force equilibrium in y-dircction requires that 
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/• 2n 

J pr r sinOdO = 2T r 


( 18 ) 


FromEqs. (17) and (18) one obtains 


T r = pr r = 


In 


(19) 


For the present paper, the total number of spokes in each rim is N s = 36 , it is seen that T r /T s = 5.73 according to Eq. (19). 
If the pretension T s in the spokes is considered as the “distributed” load, then the compressive force T r in the rim is 
equivalent to the “concentrated” load. In general, the effect of “distributed” load is much less than the “concentrated” load if 
the summation magnitude of the former is equal to the magnitude of latter. Since the effect of pretension T s is to raise the 
stiffness of spokes and that of compressive force T r is to reduce the stiffness of the rim, it is evident that the overall effect of 
increasing the pretension T s will reduce the overall stiffness of the sub-structural system composed of hub (disk), spokes and 
rim. This is the reason why the lowest several natural frequencies of the bike decrease with the increase of pretension T as 
one may see from the subsequent numerical examples. 



Fig. 3 Free-body diagram of the half rim subjected to uniform central FORCE p PER unit 

LENGTH 


2.5 Global coordinates for the two nodes of each spoke 

In the finite element analysis, preparation of input data is one of the heaviest tasks. For the bike studied in this paper, the total 
number of spokes in each wheel is N s = 36 . Thus, according to Eqs. (11a) and (lib), one must input the computer 72 pairs of 
data concerning the global coordinates of the two nodes ® and © of each spoke, (x l ,y 1 ) and ( x 2 , y 2 ), then one can obtain 
the transformation matrix | A ] of all spokes. For simplicity, the technique for determining the values of ( x t , y, ) and 
(x 2 ,y 2 ) are presented in this subsection. 

For the 36 spokes in each wheel, each of them is modeled by a two -node P-P beam element. The 1st node ® of each spoke is 
located on the rim, however, the 2nd nodes © for one half of the 36 spokes are on the 1st (hub) disk and those for the other 
half of the 36 spokes are on the 2nd (hub) disk as shown in Fig. 4 and Table 1. Because the 36 nodes are uniformly 
distributed on the rim (or hub disks), the subtended angle between any two adjacent nodes is given by A(? = 360° / N s = 10° . 
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For convenience, the 36 nodes on the rim (or hub disks) are denoted by 1, 2, .... 36 and the associated spokes by (1), (2), 
(36), respectively, beginning from 0, = 0 (or 0, = 0) as one may see from Fig. 4. It is noted that the numbering for node ® (on 
the rim) of /-th spoke is i, however, the numbering for node © (on the hub disks) of the same ;-th spoke is j as shown in 
Table 1. From Table 1, one sees that the 2nd nodes © of the spokes with “odd” numberings are on (hub) disk 1 and those 

with “even” numberings are on (hub) disk 2, and, in Fig. 4, the odd numbering spokes are denoted by the solid lines ( ) 

and the even numbering spokes by the dashed lines ( ). 


Based on the last descriptions and Fig. 4, the global coordinates for the two nodes of the arbitrary /-th spoke are given by 


with 
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(20b) 
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where ( x h , y h ) are the global coordinates of hub center H, r r and r h are the radii of the rim and the hub disks, respectively, 
while 0 X , and 0, j are the angles between the radii at nodes ® and © and the “negative” x -axis, respectively. From Eqs. 
(20a) and (20b), one sees that x 2 . = x h and y 2 : = y h if the radius of the hub disks is very small so that r h a 0 . 

Since the stiffness of the hub (disk) is much greater than that of each of the spokes, it is reasonable to assume that the hub 
(disk) is a “rigid body”, for simplicity. Furthermore, since the bike wheel is to rotate about its central axle, the center of the 
hub must be pinned. Based on the last assumptions, the translational displacements of node © of any /-th spoke are identical 
to those of the hub center //, i.e.. 


U 4 J = U yM 


(22a) 

(22b) 


Eqs. (22a) and (22b) mean that the “numberings” for the degrees of freedom of 2nd node © of each spoke are the same as 
those of the hub center H. However, the global coordinates (x 2j ,y 2j ) given by Eq. (20b) required for the determination of 
the transformation matrix [A], of the /-th spoke are different from each other for each of the spokes. 

From Fig. 4 one sees that the “moment arm” for the tensile force in each spoke with respect to the hub center H increases 
with the increase of radius r h of the hub disk, so does the restoring moment induced by the spoke tension. Thus, the effective 
stiffness of the sub-structural system composed of the rim, spokes and hub (disks) increases with the increase of r h , this is 
the reason why the lowest several natural frequencies of the entire bike structure increase with the increase of r h as one may 
see from the numerical examples given in the latter section 5. 
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Table 1 

The numberings for nodes ® and © of each spoke i ( / = 1 — 36 ) in a rim. 


Spokes (i) 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

Node ® ( i ) 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

Node © 

*Odd 

13 


27 


17 


31 


21 


35 


25 


3 


29 


(j) 

*Even 


26 


16 


30 


20 


34 


24 


2 


28 


6 


Spokes ( i ) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

Node © (i) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

Node © 

*Odd 

7 


33 


11 


1 


15 


5 


19 


9 


23 


(j) 

*Even 


32 


10 


36 


14 


4 


18 


8 


22 


12 


*Note: The 2nd nodes © for the spokes with “ odd ” numberings are on hub disk 1 and those with “even” 
numberings are on hub disk 2 (cf. Fig. 4). 



Fig.4 ( a ) Side view and ( h ) front view of the bike wheel “with hub pinned”. The coordinates of 

THE 2ND NODE © OF THE i TH SPOKE ON THE RIGID HUB DISK ARE GIVEN BY x 2 , = x h - r h COS0 2 . AND 
y 2 i = y h + r h sin 0 2 , WITH 0 2 i DENOTING THE ANGLE BETWEEN THE RADIUS r h OF THE HUB DISK// AT 

NODE © AND THE “NEGATIVE” X -AXIS. 
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III. Natural frequencies and mode shapes of the bike 

Based on the stiffness matrix [k] q , mass matrix [ m ] q , geometric stiffness matrix | k G \ tj and the transformation matrix \X\ q , 
one may obtained the property matrices of various beam elements with respect to the global coordinates and assembly of the 
latter will define the equation of motion of the entire bike structure 


[M] NxN {u} N x1 +[^]^{«}^ x1 =0 

After imposing the boundary conditions with all the constrained dof’s eliminated, one obtains 

[M]~ ~ { zT } ~ + [/T]~ ~{m}~ =0 

L J NxN l J Nx 1 L J NxN 1 J Nxl 


(23a) 


(23b) 


In the above two equations, [M] (or [ M ] ) is the overall mass matrix, [ K ] (or \K \ ) is the overall stiffness matrix, | u] (or 
[m ] ) is the overall displacement vector and [ ;T] (or [it] ) is the overall acceleration vector. Furthermore, N and N are the 
total dof’s (including unconstrained and constrained ones) and the total unconstrained dof’s, respectively. 


For free vibrations, one has 


{u(t)} = {U}e i 


(24) 


where {U } is the amplitude of {u (t ) } , CO (rad/sec) is the natural frequency of the vibrating system and / = -J — 1 . 
Substituting Eq. (24) into Eq. (23a) leads to 

([K]-cd 2 [M]){U} = 0 
Similarly, from Eq. (23b) one may obtain 

([K]-co 2 [M]){U} = 0 


(25a) 


(25b) 


Eq. (25a) or (25b) is a standard characteristic equation, from which one may determine the r-th natural frequency C0 r and the 
corresponding mode shape {C/ },. (or {I/},.), r = 1,2,3,... , of the vibrating system. In this paper, the Jacobi method [13] is 
used to solve Eq. (25), and the rth mode shape is determined by using the following global coordinates of all nodes 

X^=x i+ Ul r) (i = \-n ) 


Y^ = y i+ U^ (i = l-n n ) 


(26a) 

(26b) 


where ( x, , y, ) are global coordinates of node i required for the determination of transformation matrix | A \ q of the 
associated beam element as one may see from Eq. (11), and U 1 '.’ are displacement components of node i in the X 
and y directions, respectively, and n n is the total number of nodes of the entire vibrating system. It is noted that U {' ' and 
U 1 ' 1 are parts of the components of the r-th mode shape [ U } r determined from Eq. (25b), and the displacement components 
of the “constrained” dof’s are equal to zero. In this paper, the original configuration of the bike is obtained from the global 
coordinates ( x ; , y i ) with i = l — n n and the r-th mode shape is obtained from the global coordinates ( X* r) , Yi (r> ) with 
* = !-«,• 
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IV. Reactive forces and member forces 


For a static “constrained” structural system, its force-displacement relationship is given by 

[£]{«} = {F} 


( 27 ) 


where [ K ] and {it } have been defined in Eq. (23b), and { F } is the associated external force vector. Once the external load 
{F} is given, then from Eq. (27) one may obtain the node displacement vector 


{t7} = [£r'{F} 


(28) 


If the ;-th dof of the structural system is constrained, then the associated reactive force (or moment) is determined by 


A = I>" 

7=1 


u 

y j 


(29) 


where k /f are the matrix coefficients of the unconstrained stiffness matrix [F] - x - defined by Eq. (23a), ft , are components 
of the displacement vector {iJ} — xJ composed of the non-zero components [u } ~ given by Eq. (28) and the zero components 
corresponding to the “constrained” dof s. The reactive forces R j obtained from Eq. (29) are with respect to the global xy 
coordinate system. 

Based on the node displacement vector {w}- xl , one may obtained the node force vector {F} (l) of the ,v-th member (or 
beam element) with respect to the global xy coordinate system 

{F} W =[fc] W {M} W (30) 

where [k ] <vl and {«}' '* are the stiffness matrix and displacement vector of the ,v-th beam element. The former \k | ' ) has 
been determined before it is used to assemble the overall stiffness matrix | K \ N/N and the latter {i<} (l) is part of the 
components of {u }- . In structural design, one requires the force vector {F}"' of the ,v-th member with respect to the local 
xy coordinate system. In such a case, the values of {F} 1 ' 1 may be obtained from 

{F} (s) =[/t] (s) {F} u> 


(31) 


where [F]' ) is the transformation matrix of the ,v-th beam element. 


V. Numerical Results and Discussions 

The finite element model studied in this paper is shown in Fig. 5. For convenience, the entire bike is subdivided into three 
subsystems: the bike body, the front wheel and the rear wheel. The bike body is composed of 40 beam elements connected by 
35 nodes, either front wheel or rear wheel is composed of 36 spokes, 36 rim elements, one tire and one hub. The diameter of 
each spoke is 2 mm and is modeled by a two-node pinned-pinned (P-P) beam element. The approximate cross-section of 
each rim is to take the form as shown in Fig. 6, its sectional area is A r =60.61 mm 2 , moment of inertia of A r about its 
neutral axis n-n is I r = 1362.3259 mm 4 , the average radius of rim (based on the neutral axis n-ri) is r r = 277.5 mm (cf. Fig. 
6), each rim element is modeled by a clamped -clamped (C-C) beam element. The mass of each tire is assumed to uniformly 
distribute along the circumference of the attached rim and combined with the rim so that the effective mass density of the rim 
is given by p r = p rim + p mc , where p lim = 7850 kg/m 3 is the mass density of the rim material itself and p urc is the equivalent 
mass density of the tire with respect to the rim volume and given by p nrc = m, /V nm with m , denoting total mass of each tire 
and V rim denoting total volume of each rim. For the present example, one has m t « 1.0 kg and V rim « 0.10663 xlO -3 m 3 , 
thus, p Ae ~9378kg/m 3 . The spring constant for each tire is assumed to be k t =W/S = 100 x 9.8/0.002 = 4.9xl0 5 N/m, 
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where W = 100x9.8 N is the weight of a 100 kg rider and 8 = 0.002m is the deflection of the tire when it is subjected to the 
load of the rider. For convenience, the stiffness k of the spring-damper unit for front or rear suspension system is assumed to 
be k = 5.0x10 s N/m. The outer diameters of various tubes for the bike body are shown in Table 2, for convenience, the 
thickness of each tube is assumed to be 2.0 mm. The masses of some attachments are shown in Table 3. 

From Fig. 5 one sees that the total number of beam elements for the entire (structural) system is 186 and the total number of 
nodes is 109. Thus, the total degrees of freedom (dof’s) is 325 and the total unconstrained dof’s is 321 with the final 4 dof’s 
of the two nodes attaching the ground constrained. Besides, unless particularly stated, the numerical results of this paper are 
based on the assumption that Young’s modulus E = 2.068xl0 n N/m 2 , mass density of metal materials p = 7850 kg/ m ! , 
and there is no rider on the bike. 


Table 2 

T he outer diameters of various tubes for the bike body each with thickness 2.0 mm 


Name of tubes 

Outer diameters (mm) 

Head tube 

22 

Headset bearing tube, front fork tubes 

34 

Top tube 

28.5 

Down tube 

34 

Seat tube 

28.5 

Seat post, seat stay tubes, chain stay tubes 

25.4 


Table 3 

The masses of some parts 


Item 

Name of parts 

Mass (kg) 

Location 
(Node No.) 

Remarks 

1 

Handlebar 

0.4 

1 

*Mass of rider isn’t considered, yet. 

2 

Saddle 

0.5 

18 

*Mass of rider isn’t considered, yet. 

3 

“Chain 

0.3/2 = 0.15 

23, 34 

Uniformly shared by Nodes 23 and 34 

4 

Front hub and axle 

0.7758 

35 


5 

Rear hub, axle and 
attached sprockets 

1.5758 

34 

Combined with item 3 for node 34 

6 

Crank set, pedals and 
attachments 

1.5 

23 

Combined with item 3 for node 23 

7 

b Two tires 

1.0x2 

— 

Distributed on the associated rims 


“Summation 

7.0516 

— 

— 


“The chain mass is assumed to be equally shared by the rear axle and the crank shaft. 
b The mass of each tire is considered as the distributed mass on the associated rim. 
c Some of the attachments (e.g., the braking system ) are not included. 
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Fig. 5 The finite element model of the bike (with r h = 0 ) studied (cf. Fig. 4 for the case with r h * 0 ) 
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Fig. 6 The approximate cross-section of rim (uniform thickness: 1.1 mm ). 


5.1 Influence of attachments 

The influence of lumped masses of attachments (as shown in Table 3) on the lowest five natural frequencies of the bike 
without rider is shown in Table 4. From the table one sees that the lumped masses of attachments (no matter whether masses 
of the two tires being included or not) significantly affect the lowest five natural frequencies of the bike. This result is under 
our expectation, because the ratios between the lumped masses of the attachments (or the two tires) and the total mass of the 
186 beam elements are high. From Table 3 one sees that the total mass of the attachments (excluding the two tires) is 5.0516 
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kg, that of the two tires is 2.0 kg, and from the computer output one sees that the total mass of the 186 beam elements is 9.98 
kg ( «10 kg). It is seen that the mass ratios of the above-mentioned attachments (or the two tires) to the structural members 
(contributed to the stiffness of the entire bike) are as high as about 50% (or 20%). For this reason, in the subsequent studies 
in this paper, the effect of all attachments (including the two tires) as shown in Table 3 are taken into consideration. It is 
noted that some of the attachments (e.g., the braking system) are not included in Table 3. 

Table 4 

Influence of lumped masses of attachments on the lowest five natural frequencies of the 

BIKE WITHOUT RIDER 


Conditions 

Natural frequencies, O) (rad/sec) 

co \ 

a> 2 

<i> 3 

OJ 4 

0) 5 

No attachments 

223.37730 

260.44416 

404.32388 

455.83124 

798.58902 

With attachments 
(including two tires) 

168.05399 

207.60573 

274.94635 

311.38074 

587.16357 

With attachments 
(excluding two tires) 

178.76463 

216.50521 

302.74445 

383.54492 

647.78378 


5.2 Influence of radius of hub disk 

In most of the existing bikes, the radius of hub disk is r h ~ 25.0 mm, therefore, all the numerical results of this paper are 
based on the last value of r h . In order to study the influence of r h on the free vibration characteristics of the bike, four 
values of r h are investigated in this subsection: r h = 25.0, 15.0, 5.0 and 0.0 mm. The lowest five associated natural 
frequencies are shown in Table 5. From the table one sees that decreasing the radius of hub disk ( r h ) will reduce the lowest 
several natural frequencies of the bike as has been shown in previous subsection 2.5 and the influence is most significant on 
the 1st natural frequency ( fit, ) and decreases with increasing the vibration modes. It is seen that the influence of radius of 
hub disk ( r h ) on the 5th natural frequency ( a> 5 ) is very small so that the corresponding 5th mode shapes for the four values 
of r h look similar (cf. Figs. 7( a s ) and ( b s )). 

To save the space, only the lowest five mode shapes for r h — 25.0 and 0.0 mm are shown in Figs. 7( a t — a 5 ) and ( b x — b 5 ), 
respectively. In which. Figs. l(a l -a 5 ) denote the lowest five mode shapes for r h =25.0 mm, and Figs. l(b i -b i ) denote 
those for r h = 0.0 mm. From Fig. 7( a ) one sees that the 1st mode shape is major in (vertical) up-and-down vibration of the 
sub-structure in the “rear” suspension system (composed of the seat stays and chain stays) with respect to node 23. This is a 
reasonable result, because node 23 is a pivot and there exist a linear spring (cf. beam element No. 26 in Fig. 5) with its 
stiffness (k = 5.0x10 s N/m) much smaller than the stiffness of the other beam elements in the “rear” suspension system. Figs. 
7(a 2 ) and ( a , ) reveal that the 2nd and 3rd mode shapes are major in the heave and pitch motions of the entire bike, 
respectively. Fig. 7( a 4 ) reveals that the 4th mode shape is major in (vertical) up-and-down vibration of the sub-structure in 
the “front” suspension system, this is because there exists a soft linear spring (cf. beam element No. 5 in Fig. 5) in the front 
fork. The 5th mode shape shown in Fig. 7( a 5 ) is more complicated, because the vibrations of both the front and rear 
suspension systems are coupled. It is noted that the 1st natural frequency for r h = 0.0 mm is near zero (CO x =0.00037 
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rad/sec) as one may see from Table 5, this leads to the corresponding mode shape shown in Fig. 7( b ] ) looking like a rigid- 
body motion in (horizontal) front-and-back direction. 


Table 5 

Influence of radius of hub disk, r h , on the lowest five natural frequencies of the bike 


Radius of hub disk 
r h (mm) 

Natural frequencies, (O v (rad/sec) 

CO, 


co 3 

co 4 

0) 5 

25.0 (cf. Table 4) 

168.05399 

207.60573 

274.94635 

311.38074 

587.16357 

15.0 

105.91116 

175.50561 

236.21839 

282.04038 

585.10492 

5.0 

36.63954 

123.17002 

226.12836 

278.43011 

584.15505 

0.0 

0.00037 

114.78185 

225.45303 

278.07291 

583.95400 



Rh=25.0mm 


Rh=25.0mm 

Original configuration 


Original configuration 

3rd mode shape 


4th mode shape 
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Rh=0.0 mm 

Original configuration 

5th mode shape 



Fig. 7 (b x - b 5 ) (Continued) 


5.3 Influence of pretension in spokes 

A bicycle wheel is composed of rim, hub, spokes and tire. Among the components, the rim and hub are connected together to 
constitute a strong structure system by relying on the pretension of the spokes. However, for simplicity, the pretension T s in 
each spoke is assumed to be zero in the conventional finite element analysis and all numerical results of this paper are also 
obtained under this assumption except those presented in this subsection. Table 6 shows the influence of pretension T s in 
each of the spokes on the lowest five natural frequencies of the bike beginning from 7 = 0 with increment A7’ t =100 N. 
From the table one sees that increasing the pretension T s will reduce the lowest five natural frequencies of the bike as has 
been shown in the previous subsection 2.4. Since a strong wheel relies on enough pretension T t in each spoke and too high 
value of T s will lead to buckling of the rim, determination of the appropriate pretension T s in each spoke should be an 
optimum problem. 


Table 6 

Influence of pretension T s in each of the spokes on the lowest five natural frequencies of 

THE BIKE 


Pretension in each 
spoke, T s (N) 

Natural frequencies, co v (rad/sec) 

0) l 

a> 2 

rt» 3 

co A 

CO, 

0 (cf. Table 4) 

168.05399 

207.60573 

274.94635 

311.38074 

587.16357 

o 

o 

166.28904 

207.23125 

274.11550 

309.22295 

587.07923 

200 

164.48525 

206.83240 

273.22357 

307.11718 

586.99563 

300 

162.64151 

206.40621 

272.26276 

305.07272 

586.91278 

400 

160.75661 

205.94937 

271.22524 

303.09923 

586.83066 
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5.4 Influence of rider mass (riding gesture) on the free vibration characteristics 

In the previous subsections, the effect of rider mass is neglected and it will be studied in this subsection. It is well known that 
the mass of a rider is shared by the saddle and handlebar depending on the gesture of the rider. In this subsection, it is 
assumed that the total mass of the rider is 75 kg and it is shared by saddle and handlebar in three riding gestures: (i) m l = 25 
kg and m 18 =50 kg; (ii) m x — 18.75 kg and m 18 =56.25 kg; (iii) m l — 12.5 kg and m 18 =62.5 kg. Where the subscripts of 
m, 1 and 18, refer to nodes 1 and 18 for the handlebar and saddle, respectively, as one may see from Fig. 5. It is noted that, 
for the last three riding gestures, the handlebar shares 1/3, 1/4 and 1/6 of the total rider mass, respectively, and the remaining 
2/3, 3/4 and 5/6 of the total rider mass are shared by the saddle, respectively. For convenience, the above-mentioned three 
cases are called cases 1, 2 and 3 in Table 7, respectively, and the situation of no rider is called case 0. From Table 7 one sees 
that, in either riding gesture, the rider mass significantly affects the lowest five natural frequencies of the bike. To save the 
space, only the lowest five mode shapes of the bike for case 1 are shown in Fig. 8. Comparing with the lowest five mode 
shapes of the bike for case 0 shown in Fig. 7( a — a 5 ), one sees that the rider mass also significantly affects the lowest five 
mode shapes of the bike. From the foregoing analyses one sees that the influence of rider mass (and riding gesture) on the 
free vibration characteristics of a bike is a complicated problem and use of the information presented in the existing literature 
with rider mass neglected should be careful. 

Table 7 

Influence of rider mass (75 kg ) on the lowest five natural frequencies of the bike 


Case 

Conditions 

Natural frequencies, CO v (rad/sec) 

(O x 

co 2 

<y 3 

co 4 

CO 5 

0 

No rider (cf. Table 4) 

168.05399 

207.60573 

274.94635 

311.38074 

587.16357 

1 

m l = 25 kg , m 18 = 50 kg 

46.94273 

72.29379 

142.25082 

182.86937 

300.60258 

2 

m l = 18.75 kg , m 18 = 56.25 kg 

45.86708 

73.52113 

151.45699 

193.60493 

300.66467 

3 

m 1 = 12.5 kg , »7 18 = 62.5 kg 

44.76827 

74.62326 

164.76774 

215.10436 

300.89367 
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(e) 


Fig. 8 Influence of rider mass (case 1) on the lowest five mode shapes of the bike 

5.5 Influence of rider gesture on node displacements and internal forces of members 

In the last subsection, it has been shown that different gesture of the rider will lead to different distribution of his mass on the 
handlebar and saddle. Of course, the gravitational forces on the bike induced by the rider are also dependent on the rider 


Page | 79 
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gesture. For convenience, in this subsection, the rider gestures are assumed to be the same as those shown in Table 7, i.e.: (i) 
F l =-25x9.8 = -245 N and F 18 = -50x9.8 = -490 N; (ii) F l = -18.75x9.8 = -183.75 N and 
F ]ti = -56.25x9.8 = -551. 25N; (iii)Tj =-12.5x9.8= — 122.5 N and F lg =-62.5x9.8= -612.5N. Where the negative sign 
(-) indicates that the force F 1 (or F 18 ) is downward. The displacements of nodes 2, 11, 23, 34 and 35 (cf. Fig. 5) for the 
above-mentioned three cases are shown in Table 8(a) and the internal forces (and bending moments) of the five beam 
elements, 3, 9, 13, 37 and 40 are shown in Table 8(b). It is noted that the node displacements ( u and u ) are with respect to 
the global xy coordinate system with positive u x agreeing with positive x -axis and positive u agreeing with positive y - 
axis (cf. Fig. 5). However, in the structural design, one requires the internal forces (and bending moments), F x , F and M , , 
at the two ends (with their node numberings ® and ©) of each beam element to be with respect to the local xy coordinate 
system (instead of the global xy system) as shown in Table 8(b). 

Based on Figs. 4 and 5, and case 1 of Table 8(b), one may obtain the free-body diagrams of beam element Nos. 3, 9, 13, 37 
and 40 in the local xy coordinate systems as shown in Figs. 9(a)-(e), respectively. From Fig. 9(a) one sees that beam element 
No. 3 is subjected a “compressive” force with magnitude F x = 31 1.9N and from Fig. 9(d) one sees that beam element No. 37 
is subjected a “tensile” force with magnitude F x = 436.01 N. It is believed that the “tensile” force in chain stays (436.01 N) 
as shown in Fig. 9(d) being much greater than the “compressive” force in the front fork (311.9 N) as shown in Fig. 9(a) 
should be in agreement with the actual situations. Furthermore, the conditions of equilibrium, y F v = 0 and 

=0, are satisfied for each of the free-body diagrams shown in Figs. 9(a)-(e). In which, the length l i for the ;-th 
beam element is determined by £. = yj(x 2 -x u ) 2 + (y 2 , -y,,) 2 (i = 3, 9, 13, 37 or 40), with (x u ,y u ) and (x 2i ,y 2i ) 
representing the global coordinates of its node ® and node © as shown in Fig. 2, respectively. It is noted that, in 

Table 8(a) 

Influence of rider’s gesture on the node displacements 


Cases 

Displacements 


Numberings of nodes 



(m) 

2 

11 

23 

34 

35 

1 

F l = -245 N 
F lg = -490 N 

u x 

-0.2854E-03 

-0.2470E-03 

-.3349E-03 

-.3380E-03 

.3180E-03 

u y 

-0. 1626E-02 

-0.1683E-02 

-.1600E-02 

-.9467E-03 

-.6238E-03 

2 

F, = -183.75 N 
F lg = -551 .25 N 

«, 

-.3909E-03 

-.2914E-03 

-.3031E-03 

-.3064E-03 

.2865E-03 

Uy 

-.1539E-02 

-.1708E-02 

-.1610E-02 

-.1005E-02 

-.5656E-03 

3 

Fj = -112.5 N 
F 18 = -612.5 N 

u x 

-.4964E-03 

-.3358E-03 

-.2713E-03 

-.2749E-03 

.2550E-03 

Uy 

-.1452E-02 

-.1734E-02 

-.1619E-02 

-.1063E-02 

-.5073E-03 
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Figs. 9(a)-(e), the counterclockwise (CCW) moments (M _ u or M , 2i ) are positive because they are the bending moments 
about the positive z (or Z )-axis. Similar to Figs. 9(a)-(e), one may obtain the free-body diagrams of the beam element Nos. 3, 
9, 13, 37 and 40 for the other cases of Table 8(b). 

The configurations of the bike after deformations for the above-mentioned three cases (cf. Table 8) are shown in Figs. 10(a)- 
(c), respectively. In each case, all deformations (or displacements) of the nodes induced by the external forces F l and F lg 
are normalized by the maximum one in that case and then multiplied by a factor 0.05 to obtain the appropriate dimensions of 
the deformed bike for plotting. Because the external loads, F t and F lg , are downward and Fi < F ls , all the vertical node 
displacements ( u ) as shown in table 8 are in negative (-) y -direction and all the horizontal node displacements ( u x ) as 
shown in same table are in negative (-) x -direction (except those for node 35). It is the last reason, the deformations of the 
entire bike for the three cases have the same trend of tilting leftward, so that the deformed configurations of the bike for the 
three cases shown in Figs. 10(a)-(c) look similar. However, because the vertical force on the handlebar ( F { ) is maximum in 
case 1 and minimum in case 3, so is the vertical deformation of node 2 near the handlebar (cf. Fig. 5) as one may see from 
Figs. 10 (a)-(c). 


Table 8 (b) 

Influence of rider’s gesture on the internal member forces and moments 


Cases 

Nodes 

Forces (N) 


Numberings of beam elements 




or 


(Numbering of 1 

st node ^Numbering of 2 nd node) 




moments 

3 

9 

13 

37 

40 



(Nm) 

(3-A4) 

(8-»9) 

(12-M3) 

(23-»31) 

(33->34) 

1 

® 

F x 

.31190E+03 

.46246E+03 

-.46448E+03 

-.43601E+03 

-.43601E+03 

F l = -245 N 


F y 

-.32154E+02 

.16723E+03 

-.60647E+02 

.13000E+01 

.13000E+01 

F w = -490 N 


M z 

-.17975E+02 

. 1 3434E+02 

.48521E+01 

*0 

-.40949E+00 


© 

F 

X 

-.31190E+03 

-.46246E+03 

.46448E+03 

.43601E+03 

.43601E+03 



Fy 

.32154E+02 

-.16723E+03 

.60647E+02 

-.13000E+01 

-.13000E+01 




.14380E+02 

.83103E+01 

-.12255E+02 

.13650E+00 

.54599E+00 

2 


F 

X 

.28249E+03. 

.46990E+03 

-.50805E+03 

-.48138E+03 

-.48138E+03 

F t = -183.75 N 


F 

y 

-.29705E+02 

.17371E+03 

-.73535E+02 

.14543E+01 

.14543E+01 

F lg =-551.25 

N 


M 

Z 

-.16606E+02 

. 1 2482E+02 

.21355E+01 

0 

-.45809E+00 


© 

F x 

-.28249E+03 

-.46990E+03 

.50805E+03 

.48138E+03 

.48138E+03 



Fy 

.29705E+02 

-.17371E+03 

.73535E+02 

-.14543E+01 

-.14543E+01 



M 

Z 

.13285E+02 

.10104E+02 

-.1 1 1 12E+02 

.15270E+00 

.61079E+00 

3 

® 

F x 

.25308E+03 

.47734E+03 

-.55162E+03 

-.52676E+03 

-.52676E+03 
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F l =-112.5 N 


F 

y 

-.27256E+02 

.18019E+03 

-.86424E+02 

.16086E+01 

.16086E+01 

F Ig = -612 .5 N 


M 

Z 

-.15237E+02 

.1 1531E+02 

-.581 19E+00 

0 

-.50669E+00 


© 

F x 

-.25308E+03 

-.47734E+03 

.55162E+03 

.52676E+03 

.52676E+03 



F y 

.27256E+02 

-.18019E+03 

.86424E+02 

-.16086E+01 

-.16086E+01 




.12189E+02 

.1 1898E+02 

-.99682E+01 

.16890E+00 

.67559E+00 


* Pivot at node 23. 
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Fig. 9 Free-body diagrams for Nos. ( a ) 3, ( b ) 9, ( c ) 13, ( d ) 37 and ( e ) 40 beam elements based on 
Figs. 4 and 5, and case 1 of Table 8( b ). All beam elements are C-C except No. 37 beam element 

IN ( d ) BEING P-C. 
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Case 3 (F,= -1 12.5 N, F„. -612.5 N) 

Original configuration 

Deformed configuration 



Fig. 10 The deformed configurations of the entire bike subjected gravitational forces of 

the rider: ( a ) case 1 ; ( b ) CASE 2; ( c ) case 3. 

VI. Conclusion 

Based on the foregoing numerical-analysis results, the following conclusions are drawn: 

1. Because the ratio of the lumped masses of attachments (including the two tires) to total mass of the structural members 
(contributed to the stiffness of the entire bike) is high (near 50%), the dynamic analysis of a bike should be conducted 
with the effect of its lumped masses of attachments considered. 

2. The “moment arm’’ for the tensile force in each spoke with respect to the hub center increases with the increase of 
radius r h of the hub disk, so does the restoring moment induced by the spoke tension. Thus, the effective stiffness of 
the sub-structural system composed of the rim, spokes and hub (disks) increases with the increase of r h and so do the 
lowest several natural frequencies of the entire bike. 

3. Because the compressive force T r in a rim induced by the pretension of all its spokes is much larger than the pretension 
7 . in each spoke, in spite of the fact that the pretension 7 t can raise the stiffness of each spoke and the compressive 
force T r can reduce the stiffness of the rim, the overall effect of increasing the pretension T s will reduce the overall 
stiffness of the sub-structural system composed of hub (disk), spokes and rim, thus, the lowest several natural 
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frequencies of the bike will decrease with the increase of pretension T s . Since a strong wheel relies on enough 
pretension T s in each spoke and too high value of T s will lead to buckling of the rim, determination of the appropriate 
pretension T : in each spoke should be an optimum problem. 

4. The free vibration characteristic of a bike is significantly affected by the mass of its rider. Thus, for a bike to 
accommodate various riders, some of its pertinent parameters (e.g., the axial stiffness of the “spring-damper units” for 
rear and front suspension systems) should be adjustable. In addition, most of information obtained from the bike 
without a rider in the existing literature, its applicability in practice seems to need further studies. 

5. For a bike subjected to the gravitational forces of a rider, its node displacements and internal forces (and bending 
moments) of the structural members are significantly affected by the riding gesture of the rider. Based on the theory 
presented and the computer program developed for this paper one may easily obtain the last information, this should be 
benefit for designing a save and comfortable bike. 
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Appendix A 

Stiffness and mass matrices of the P-P, P-C and C-C beam elements 

The stiffness matrix [ k \ pp and mass matrix [tn] pp for the P-P beam element are given by [9] 
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The stiffness matrix [k\ PC and mass matri x [ in ] pc for the P-C beam element are given by [10] 
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The stiffness matrix [ k \ (;c and mass matrix [tfl\ cc for the C-C beam element are given by [9] 
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